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SEQUENTIAL  TESTS  FOR  THE  DRIFT  OF  A  WIENER  PROCESS 
WITH  A  SMOOTH  PRIOR,  AND  THE  HEAT  EQUATION 


by  Gordon  Simons1,  Yi-Ching  Yao  and  Xizhi  Wu 
University  of  North  Carolina  in  Chapel  Hill, 

Colorado  State  University  and  University  of  California  at  Davis. 


Methods  are  described  which  permit  one  to  work  with  continuous— time 
optimal  stopping  problems,  using  the  heat  equation,  even  when  the  prior 
placed  on  the  drift  parameter  of  a  Wiener  process  is  not  normal.  The  details 
of  the  method  are  worked  out  for  Chernoffs  problem  of  testing  the  sign  of 
the  drift  parameter  when  the  prior  is  "smooth". 
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1.  Introduction.  Chernoff  (1961,  1972)  and  others  have  shown  how  the  heat  equation  and 
associated  free  boundary  problems  arise  naturally  as  mathematical  tools  for  approximating 
Bayes  sequential  procedures.  Typical  applications  are  concerned  with  normally  distributed 
observations  which  depend  on  an  unknown  mean  9  and  a  normal  prior  for  9.  The  present 
paper  has  two  objectives:  (i)  to  show,  very  generally,  that  the  heat  equation  is  still  relevant 
when  9  is  not  normally  distributed,  and  (ii)  to  investigate  the  problem  of  testing  the  sign 
of  the  normal  mean  for  smooth  priors.  A  smooth  non— normal  prior  can  be  viewed  as  a 
perturbed  normal  prior,  and  it  is  possible  to  very  precisely  describe  how  the  perturbation 
affects  the  associated  free  boundary. 

Since  the  subject  of  free  boundary  problems  is  fraught  with  technical  details,  a 
conscious  decision  has  been  made  here  to  emphasize  exposition  rather  than  rigor. 

2.  Bayes  problems  and  continuous  time  approximations.  Suppose  Xt,  X2,  •  •  •  are  potential 
observations  which  are  independent  and  N(0,<r2),  where  a2  is  positive  and  known,  and  6 
is  unknown.  Set  Sn  =  Xt  +  •  •  •  +  Xn,  n>0,  and,  for  definiteness,  a2  =  1.  We  shall  assume 
that  we  are  dealing  with  a  concrete  statistical  problem  which  is  expressed  within  a 
decision— theoretic  framework  with  a  suitable  loss  structure,  and  that  9  has  a  prior  density 
g.  The  Bayes  sequential  procedure  is  sought. 

Let  d(z,n)  denote  the  posterior  Bayes  risk  associated  with  stopping  at  time  n,  with 
Sn  =  z,  and  making  an  optimal  terminal  decision.  The  primary  task,  in  a  Bayes  sequential 
problem,  is  to  find  the  stopping  time  N,  within  {0,  1,  2,---},  which  minimizes  the 
expectation  E{d(Sv,N)}. 

We  are  only  concerned  here  with  the  commonly  used  approximation  based  on  the 
replacement  of  {Sn,  n>0}  by  a  Wiener  process  X  =  { X(t) ,  t>0}  with  drift  rate  9-.  Given 
9 ,  the  differential  dX(t)  has  mean  0-dt  and  variance  dt.  The  task  is  to  find  the 
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stopping  time  r  in  [0,®)  which  minimizes  E{d(X(r),r)}.  Such  problems  can  be  viewed 
as  Markovian,  with  states  (x,t),  and  the  search  for  optimality  can  be  viewed  as  a  search  for 
an  "optimal  continuation  set"  C  in  the  (x,t)  plane:  The  optimal  stopping  time  r  =  r(x.t) 
is  the  first  time  that  X  reaches  the  boundary  of  C  if  (x,t)  is  in  C.  And  it  is  simply  t 
otherwise. 

Let  b(x,t)  :=  E{d(X( r),r)},  represent  the  minimal  Bayes  risk  that  can  be  achieved 
by  stopping  optimally  starting  from  state  (x,t).  The  same  argument  as  that  used  by 
Chernoff  (1972,  p.  92)  can  be  used  to  show  that  it  satisfies  a  diffusion  equation: 


bxx  +  2v?(x,t)-bx  +  2bt  =  0,  (x,t)eC, 


where 


tfx,t)  :=  E{fl|X(t)=x}  = 


JZ  9  exPlx*  ~  M  d9 

JZ  exp{x<?-t^/2}  g (9)  d 7 


I 

Ik!? 


And,  with  the  possible  exception  of  certain  points  of  singularity,  the  boundary  conditions 


(3)  b  =  d,  bx  =  dx,  bt  =  dt 

hold  on  the  boundary  of  C.  Frequently,  equation  (1)  and  the  first  two  conditions  in  (3) 
jointly  determine  the  function  b  and  the  set  C.  But  these  conditions  do  not  always 
guarantee  a  unique  solution. 

We  note  in  passing  that  the  ratio  in  (2)  is  meaningful  for  t  >  0  even  when  the  prior 
g  lacks  a  finite  first  moment. 


3.  ChernofFs  negative  s  scale.  It  is  a  standard  practice  to  assume  that  the  prior  g  is 
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normal  with  some  known  mean  /rQ  and  variance  aQ2.  This  makes  the  marginal 
distribution  of  X  Gaussian  but  something  other  than  a  Wiener  process  with  drift.  It  does 
not  make  equation  (1)  particularly  attractive  to  work  with.  For  this  reason,  Chernoff 
(1972)  and  others  have  found  it  more  convenient  to  work  with  the  posterior  mean  process 
Y  =  { V(s),  0  <  s  <  crQ 2},  defined  by 


Y(s)  =  E{0| X(t)}  =  (X(t)  +  ^ /0/( t  +  a  -2), 


where  s  =  l/(t  +  <rQ'2)  is  the  posterior  variance  at  time  t.  Then  Y  is  a  Brownian 
motion  evolving  backwards  in  time,  as  measured  in  s  (hence  the  name  "negative  s  scale"). 
And  b(x,t),  when  described  as  a  function  of  (y,s),  satisfies  the  much  more  familiar  heat 
equation  uyy  =  2us,  and  boundary  conditions  comparable  to  those  in  (3).  A  frequent 
additional  advantage  is  that  Y(s)  is  statistically  more  relevant  than  X(t).  I.e.,  d(x,t)  is 
more  naturally  described  as  a  function  of  (y,s). 

If  the  prior  g  is  non-Gaussian,  the  marginal  distribution  of  X  is  non— Gaussian. 
Moreover,  the  posterior  mean  is  non-Gaussian,  and  it  no  longer  assumes  the  simple  linear 
form  shown  in  (4). 

The  loss  of  the  linear  form  is  unfortunate;  for  it  is  the  reason  Y  is  a  Brownian 
motion.  It  will  be  seen  in  the  next  section  that  it  is  possible  to  continue  to  work  with  a 
variant  of  this  linear  form.  We  will  be  content  here  with  pointing  out  that  the  posterior 
mean  is  approximately  equal  to  a  linear  form  when  t  is  large  and  the  prior  g  is 
sufficiently  smooth:  It  can  be  shown  that  (t+c)E{0j X( t) }  -  (X(t)+a)  is  a  convergent 
martingale,  for  arbitrary  constants  a  and  c,  whenever  g  is  continuously  differentiable,  in 
which  case  the  limit  equals  g'(0)/g(0)  -  a  +  c0  a.s.  A  linear  form  will  be  used,  as  a 
mathematical  convenience,  when  we  consider  the  problem  of  testing  the  sign  of  the  normal 
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mean.  The  results  so  obtained  will  need  to  be  restated  in  the  coordinate  system  of  the 
posterior  mean  before  meaningful  comparisons  can  be  made  with  Chernoff's  results  for 
normal  priors. 

4.  Transformations.  We  are  concerned  here  with  transformations  of  b  which  satisfy  the 
heat  equation  uxx  =  2ut.  We  start  with  a  transformation  which  satisfies  the  backward 
heat  equation  uxx  =  -  2ut: 

Theorem  1.  The  function 

(5)  b  (x,t)  :=  b(x,t)*0(x,t) 


satisfies  the  backward  heat  equation  in  C,  where 


iKx, t)  :=  /  “  exp{x0-t^/2}  g(0)  dfl. 


The  proof  follows  by  direct  calculations. 


Now  set 


d  (x,t)  :=  d(x,t)-  t/<x,t). 


The  functions  b  and  d  can  be  viewed  as  surrogates  for  b  and  d,  respectively: 
Instead  of  trying  to  solve  the  free  boundary  problem  represented  by  equation  (1)  with 

X 

boundary  conditions  (3),  one  may  solve  an  equivalent  problem  using  the  surrogates  b 


and  d  . 

What  we  have  done  can  be  understood  from  a  probabilistic  perspective  as  well:  The 

transformation  corresponds  to  a  change  of  probability  measures  with  ^X(t),t)  equal  to 

* 

the  Radon-Nikodym  derivative  dP/dP  for  the  (current)  sigma-field  <7{X(u):  u<t } . 


bt'lA'M  Lt'll  ft.h 


Under  the  new  probability  measure  P  ,  X  is  a  Brownian  motion,  d  and  b  are  the  proper 
analogues  of  d  and  b  respectively,  and  the  pro[x>r  analogue  of  (1)  is  the  backward  heat 
equation. 

To  get  from  the  backward  heat  equation  to  the  heat  equation,  let  y  =  (x  +  a)/(t,  + 
c)  and  s  =  1  /( t  +  e)  (t  >  —  c),  and  let  (V  lx'  the  image  of  C  in  the  (y,s)  plane,  where, 
for  now,  a  and  c  are  arbitrary  constants.  Thus,  x  =  y/s  —  a  and  t  —  1 /s  -  c.  Further. 

_i  i 

let  o  lx' t  he  standard  normal  density  and  n(v,s)  =  s  •’  0( y / s' ) . 

Theorem  2.  The  function 


b  (y.s)  m  b  (x.t)-n(y,s). 


satisfns  tht  hint  tqiitilioii  in  tin  rariabhs  (y.s)  within  (". 

Again  t  h<‘  proof  follows  by  direct  calculations.  And  again  a  probabilistic 

interpretation  can  be  given:  Assume  X  -  (X(t).  1>0}  is  standard  Brownian  motion  under 

the  probability  measure  P  .  and  c  >  0.  Then  the  process  Y  =  (Y(s)  =  (X(t)  +  a ) / ( 1  + 

*  * 

c ).  0*,s<l/c}  is  Brownian  motion  in  reverse  time  s.  under  a  probability  measure  P  .  if 
for  each  s.  the  Kadon-Nikodym  derivative  (IP  /(IP  for  the  (current)  sigma-field 
i7{Y(ii):  s< u <  1  /<* }  is  proportional  to  n ( A’ ( s ) .s ) .  ( Tin1  coefficient  of  proportionality  equals 
| c  •  n I  a .c )] " 1 . ) 


leorems  1  and  2  together  vield  the  lormula 


b  (  y.s)  1  >(  x .t )  •  /  ,'n(y— f/.s)-e  g(  0)  dfl. 


Again,  each  solution  of  the  original  free  boundary  problem  can  be  described. 


equivalently,  in  terms  of  surrogates,  b  and 


(9)  d  (y,s)  =  d(x,t)-/^  n(y-0,s)-e  a<?+C^/2  g(0)  d0. 

While  a  particular  solution  b,  so  found,  may  not  be  the  solution  sought  for  the  sequential 
Hayes  problem,  one  of  the  obtainable  solutions  in  the  (y,s)  plane  will,  necessarily, 
correspond  to  the  sequential  Bayes  problem;  there  is  a  simple  one-to-one  correspondence 
between  solutions  in  the  two  spaces. 

Notice  that  the  original  problem  has  been  converted  into  an  equivalent  problem  in 
the  context  of  Chernoff's  negative  s  scale1  described  in  Section  3.  This  has  several 
advantages.  An  obvious  advantage  is  that  the  heat  equation  has  been  widely  studied  by 
probabilists  and  statisticians.  Many  solutions  are  known,  and  these  can  be  used  to  obtain 
asymptotic  expansions  and  other  analytic  approximations  to  the  solution  of  the  free 
Ixnmdary  problem.  Another  convenience  is  that  Chernoff  and  Pctkau  (1.981,  1.9875)  haw 

x  x 

developed  numerical  algorithms  for  the1  negative  s  scale  which  find  b  and  (",  for  quite 
general  functions  d  .  A  final  potential  advantage  is  that  one  is  able  to  place  a  large 
family  of  optimal  stopping  problems  within  a  common  context;  each  prior  g  contributes  a 
member.  This  invite's  new  t\|>es  of  comparisons.  An  example  of  this  is  given  in  the  next 
sect  ion. 

It  is  now  a  quick  matter  to  completely  recover  Chernoff's  framework  when  the  prior 
g  is  normal  \(/<  ,c7  2):  Simply  set  a  =  f‘0/oQ2  and  e  =  l/cro2.  Then  (8)  and  (9)  become 

x  *  *x 

b  (y.s)  =  n(//0.cT02) •  b(x.t )  and  d  (y.s)  =  n(/c  .<ro2)  •d(x.t ).  respectively.  The  positive 
constant  "(//  ,cr  2)  contributes  nothing  of  significance  to  the  problem,  and  can  simply  be 

x  X 

dropped.  So.  effectively,  we  are  concerned  with  the  transformation  cl  (y.s)  =  d(x.t). 
where  y  =  (x  +  //  /cr  2)/(t  +  (JQ  2)  and  s  =  l/(t  +  a  2).  and  we  seek  a  suitable  solution 


/  < 


b  (y.s)  of  the  heat  equation.  Chernoff,  by  focusing  attention  on  the  posterior  mean 
process  (4),  reaches  the  same  context. 

More  generally,  suppose,  for  now,  the  prior  g  is  strictly  positive  and  smooth  in  the 
sense  that  logg(0)  is  expressible  in  a  Taylor  series  expansion  : 

(10)  log  g(9)  =  ao  +  ai#  +  a202/2  +  a.393/6  +  •  •  •. 

It  is  convenient  to  think  of  aj  and  —  a2  as  prior  parameters  with  all  the  other  subscripted 
a's  held  fixed  except  ao,  which  must  be  chosen  to  make  the  integral  of  g  equal  to  one.  Then 
the  posterior  density  of  g  given  X(t)  =  x  remains  in  the  family,  and  it  has  posterior 
parameters  x  +  ai  and  t  -  a2.  This  is  a  generalization  of  the  normal  family,  which 
corresponds  to  ai  =  oQ2  and  a2  =  -  1  /<ro2,  with  ak  =  0  for  k  >  3.  In  general,  the 
integrals  in  (8)  and  (9)  take  the  form 

(11)  exp(ao)-  /_^  n(y-0,s)  -exp  {a3^/6  +  a4^/24  +  •  •  • }  d# 

when  one  sets  a  =  ai  and  c  =  —  ai-  This  makes  y  =  (x  +  a()/(t  —  a2). 

5.  Testing  the  sign  of  the  drift. 

** 

5.1.  The  function  d  (y,s).  Assume  the  cost  of  sampling  is  one  per  unit  time  and  the  cost 
due  to  a  wrong  assessment  of  the  sign  of  the  drift  rate  9  is  2k|0|,  where  k  >  0.  This 
leads  to 


d(x.t)  =  t  +  2k -min 


/  10  0ex#~t<B/2g(0)  dd  / 


0 

— oo 
oo 

— X 


|,|ex*-l*/2g(,)  dd 

ex0-t^/2  g(0)d0 


d**(y.s)  =  (s'1  +  a 2)/^  n(y-0,s)e~ai*-a2^/2g(0)  d 9  + 

2k ■  min j/“  0-n(y-0,s)e_ai^“a2^/2g(0)  d 9-J_^  9-n(y-9,s)e~^$~^ !2g{9)  d9 


Before  proceeding,  we  shall  find  it  convenient  to  subtract  away  a  solution  of  the  heat 

equation,  /  ”  (a2+k|  0|  )n(y— 0,s)e  a'^“a2^/“g(^)  d0,  and  then  tc  'fivide  by  the  positive 
constant  e30,  which  will  not  affect  the  boundary  of  C’.  (This  particular  solution  is  defined 
for  0  <  s  <  oo  if  a2  >  0,  and  for  0  <  s  <  — a2~*  if  &2  <  0.)  Keeping  the  same  notation,  as 

a  convenience,  we  obtain 


d**(y.s)  =  s-‘/4  n(y-0.s)e_a<)~ai*“a2^/2g(0)  d9 

(12) 

-k-|/_£  9-n(y-9.s)e~^l9~^/2g(9)  d9\. 

Using  the  fact  that  the  function  t!’x(x,t)  =  /_^  0-exp{x0  -t^/2}  g( 9)  d9.  is 
monotone  in  x,  one  can  show  that  the  second  integral  in  (12)  vanishes  along  a  line 
(y  (s).s).  It  is  positive  above  the  line  and  negative  below.  It  can  be  argued,  along  the  lines 
of  Sobel  (1953),  that  the  continuation  region  C'  assumes  the  form  {(y.s):  y'(s)  <  y  < 
y*(s).  s  >  0}  with  (vo(s),s)  inside  C'.  One  of  our  main  objectives  is  to  describe  the 
asymptotic  form  of  the  boundaries  y±(s)  as  s  goes  to  zero.  It  can  be  shown  that 


alp)  =  -  a3^3/2  -  a,5/95/8  -  (a?  +  20a3a4)p7/48  -  0(p9)  as  s  -  0, 


AA  A  -•  A  ■-  -  ,>.V,V.AV.V,  AA,  AA, -V, 


I 


where  aQ(p)  :=  y0(s)s  2  and  p:=s!.  This  function  is  identically  zero  when  e  l  g{6) 
is  an  even  function,  and,  hence,  for  all  normal  priors. 

**  _ i 

5.2  Expansions  of  d  .  Let  a  :=  ys  2  and  let  r\p  {a)  be  the  r-th  moment  of  a  standard 

normal  random  variable  about  -  a,  i.e.,  pr(a)  =  /_^  0(u)(u+a)rdu/r!.  Thus  po(a)  =  1, 

Pi(a)  =  and  pT(a)  =  +  a-pT.x{a))/ r  for  r  =  1,  2,  •••.  It  is  useful  to  know 

that  p T{ot)pv  satisfies  the  heat  equation  in  the  variables  (y,s). 

*  * 

A  preliminary  step  is  to  expand  the  function  d  in  terms  of  thn  pt(ct)'s  and 
powers  of  p:  For  bounded  o,  as  p  goes  to  zero. 


f_£  n(y— 0,s)-e“^~ai^a2^/2-g(0)  dO 

~  /  A  n(y— 0,s)-e_a<3_ai^a2^/2.g(0)  dd 
\9\<P' 

=  /  A  n(y-0,s ).ea3^/6+a4^/24+-"  d0 

\e\  <p' 

=  J  ,  n( v— 0,s)  •  { 1  +  a3^/6  +  a4^/24  +  •  •  •}  dO 
\0\<P' 

~  /_£  n(y~M  ■  { 1  +  a3^/6  +  a4^/24  +  •  •  • }  dO 
=  c0/io(a)  +  Ci/ii(a)/>  +  c2/i2(Q')p2  +  C3/i3(o  )pP  +  ■  •  •, 


where 


Co  =  1,  Ci  =  c2  =  0,  c3  =  a3,  c4  =  a4,  c5  =  a5,  C6  =  a6  +  10a32. 
C7  =  a7  +  35a3a4,  cs  =  as  +  35a4“  +  56a3a5, 

Cg  —  ag  +84a3a^  ■+■  126a4a5  +  280a32, 

C10  =  ato  ■+■  120a3a7  +  126a5_  +  210a4a6  +  2100a3-a4. 


The  deletion  of  the  set  {|0|  >  p2}  from  the  range  of  integration,  and  its  subsequent 
reinsertion,  contribute  negligible  terms  of  order  e  Likewise,  one  obtains 


P 

& 


f_£0-n(y-0,s)-e  *°  (0)  dO 

~  c oM<*)p  +  2ciii2(a)p2  +  3 c2/i3(a)^3  + 


and,  hence, 

(13) 


d  (y.s)  -  ^0Mp'2  +  s  (wW®)  *  rkcr.,/ir(a)]/?r, 
r>l 


where  the  sum  includes  a  finite  number  of  consecutive  terms.  The  upper  and  lower  signs  in 

"  r  "  are  applicable  above  and  below  the  line  aQ{p),  respectively. 

Note  that  the  even  and  odd  functions  of  a  always  go  with  the  even  and  odd  powers 

** 

of  /?,  respectively.  The  same  must  be  true  for  the  expansion  of  b  .  The  presence  in  (13) 
of  a  negative  second  power  of  p  forces  us  to  say  something  more  about  solutions  of  the 
heat  equation. 

5.3  Special  solutions  of  the  heat  equation:  For  any  integer  r,  V T{a)pr  is  a  separable 
solution  of  the  heat  equation  (in  y  and  s)  if  Vr  satisfies 

(14)  V”(a)  +  a-V;(a)  =  r-Vr(a). 

The  even  and  odd  solutions  of  (14)  are,  respectively: 

G,(  a)  =  1  +  \ -o?  +  ^  „<  +  •  •  ■  +  r(t-2) a* 


Hr(a)  =  o  +  +  (r  MU  3)-a5+  ■■■  +  tlTHl 


— 1  )(r— 3 )  •  •  -(r— 2n+l : 


>n+l  ! 


a2"*1  + 


The  function  pr(a)  is  a  multiple  of  Gr(a)  when  r  is  an  even  integer,  and  a 
multiple  of  Hr(a)  when  r  is  odd  (r  =  0.1.2. •  ■  •)•  For  notational  convenience,  we  will  let 
/i.2(o)  =  G.j(a)  =  1  -  2o2/2  +  2-4q4/24  -  •  •  ■  and  /m(o)  =  H.,(o)  =  a  -  2a3/6  - 
2-4os/120 - . 

mm 

We  shall  seek  an  asymptotic  solution  for  b  of  the  general  form 


b  (y.s)  -  E  ur/ir(a)pr, 
r>— 2 


where,  again,  the  sum  includes  a  finite  number  of  consecutive  terms,  and  the  ur's  are 


constants  which  are  to  be  determined  from  the  boundary  conditions. 


5.4  Four  boundary  conditions:  Two  boundary  conditions  are  described  by  equating  the 

mm  ** 

expansion  of  d  with  the  presumed  form  of  the  expansion  of  b  .  This  yields 


(u.j/i-jfo)  -  Dp'2  +  u ,xp.x(a)p''  +  u^fo) 


+  E  [( ur  ±  rkcrM)u  (a)  -  c n2pT.2(a)]pT  =  0, 
r>  1 


where  "  *  "  refers  to  the  upper  and  lower  boundaries  y±(s),  respectively.  This  is  to  hold 

■  I  _ i  **  »* 

when  a  =  a  (s)  :=  y  (s)s  ’.  The  further  requirement  that  dv  and  by  should  be 
equal  on  the  boundaries  can  be  incorporated  by  differentiating  (16)  with  respect  to  a: 


u-2 p-2'{a)p'2  +  u-i p-i(a)p~l  +  £  [(ur  ±  rkcr.,)#ir'(o)  -  cr+2/ir+2'(a)]/)r  =  0. 

r>  1 


Since  pr'  =  /ir.,  (r  =  -1,  1,  2.  3,  •  •  • ),  this  becomes 


S5WS»S^S 


-  13  - 


lj  L|*W| 


j  wax* 


a»  «■*  A  «>a  **1.'^ 


(IT)  u-2/i-3( a)^'2  +  u-i/MtaJ/T1  +  S  [(ur  ±  rkcr.,)/zr.,(Q) -cr.2/irt,(a)]pr  =  0, 

r>  1 


where  p-z(a)  ■=  p-2 '{a)  =  —2a  +  2-4o3/6  —  2-4,6o5/120  +  •  •  •• 

5.5  The  expansions  of  a±.  Following  Chernoff  (1972,  p.  99),  we  anticipate  that  the 
boundary  functions  q±(s)  have  a  general  asymptotic  form 

( 18)  a±(s)  -  b3±p3  +  b4±p4  +  b5 V  +  •  •  • . 

The  constants  bj^  and  the  coefficients  ur  are  determined  by  the  boundary  conditions. 
The  details  are  messy  but  the  idea  is  quite  simple:  One  replaces  a  everywhere  in  (16)  and 
(17)  by  the  presumed  form  of  the  asymptotic  expansions  of  a^s)  given  in  (18).  This 
yields  four  series  in  p  which  can  be  viewed  as  being  identically  equal  to  zero.  By  setting 
the  coefficients  of  the  powers  of  p  equal  to  zero,  one  obtains  a  system  of  equations  which 
can  be  solved  to  yield  the  coefficients  ur  (r  =  -2, -1,0,*  •  •)  and  bj  (j  =  3,4,-  •  •).  One 
finds: 

U-2  =  1,  u-i  =  uo  =  0,  ui  =  -  aa/2,  U2  =  a4/4,  U3  =  -  as/4, 

u4  =  (ag  —  2a32  -  12k2)/6,  U5  =  -  (a7  +  lla3a4)/6,  U6  =  (ag  +  8a3as  +  35a42)/8. 

u7  =  -  (ag  +  36a3ag  +  72a33  +  78a4a5)/8  -  38k2a3, 

u8  =  (aio  +  40a3a7  +  66as2  +  210a4a«  +  380a32a4)/10  -  108k2a4, 

b4±  =  b6±  =  b8±  =  •••=  0,  b3±  =  -  a3/2  ±  k/2, 
b5±  =  -  a5/8,  b7±  =  -  (a7  +  20a3a4)/48  ±  ka4/4, 

b9±  =  -  (a9  +  56a3a<5  +  96a4a5  +  96a33)/384  ±  k(a«  +  4a32)/l6  +  k2a3/24  r  k3/6. 


b 


5.6  Expansions  in  the  coordinate  system  of  the  posterior  mean.  In  order  to  clearly 
interpret  (and  simplify)  the  results  described  in  the  previous  subsection,  it  is  necessary  to 
change  from  the  coordinate  system  based  on  the  mapping  y  =  (x  +  ai)/(t  —  a2)  to  the 
coordinate  system  based  on  y  =  E(0jX(t)=x)  =  ipx(x,t)/ip(x,t).  It  is  easily  checked  that 
y  is  a  continuous  and  strictly  increasing  function  of  x  when  t  is  held  constant.  So,  there  is 
a  well-defined  one-to-one  mapping  from  (y,s)  to  (y,s).  One  obvious  simplification  is  that 
the  function  yQ(s)  gets  mapped  into  the  the  s-axis  of  the  (y,s)  coordinate  system.  The 
most  interesting  question  is:  "How  are  y±(s)  mapped?"  This  is  easily  answered  using  the 
information  at  hand:  Let  y  (s)  denote  the  mappings  of  y  (s).  The  asymptotic 
expansions  of  y±(s)s~’  as  p  goes  to  zero  take  the  form  ba^/?3  +  bs^/?5  +  b p1  +  bg±^9 
+  •  •  • ,  where 

ba^  =  ±  k/2,  bg^  =  0,  bT1^  =  ±  ka4/2,  b^  =  ±  k(a6  -f  4a32)/8  +  k2aa/6  *  k3/6. 
Several  observations  can  be  made: 

(i)  Chernoff's  expansion  for  a  normal  prior  has  the  same  first  term  (also  obtained  by 
Bather  (1962)).  Our  second  nonzero  term  does  not  show  up  in  his  expansion  since 
normal  priors  have  a4  equal  to  zero.  His  second  nonzero  term  is  b^  =  *  k3/6. 
(Note  that  our  k  is  one— half  of  Chernoff's  k.) 

(ii)  The  evidence  for  a  non-symmetric  optimal  stopping  boundary,  in  the  (y,s)  coordinate 

system,  first  appears  in  the  formula  for  bg  ,  and  then  only  if  a3  is  not  zero.  Its 
effect  is  to  shift  the  continuation  region  upward  when  positive,  and  in  the  opposite 
direction  when  negative. 

(iii)  The  coefficients  a4  and  a6  produce  similar  effects:  When  positive,  each  causes  the 

continuation  region  to  be  enlarged. 


mili'tiii  iCiU' 
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(iv)  It  seems  somewhat  surprising  that  a4  shows  up  earlier  than  a3  in  the  expansions  of 
y±,  and  a6  earlier  than  as  (which  presumably  appears  in  bit).  It  seems  that 
departures  from  normality  more  readily  affect  the  size  of  the  continuation  region 
than  the  location  of  its  midline,  at  least  when  s  is  small  (large  t). 

We  are  uncertain  how  to  intuit  observations  (ii) — (iv). 

It  should  be  pointed  out  that  we  have  assumed  much  more  about  the  prior  density 
g  than  is  necessary.  Clearly,  for  the  results  we  have  described,  one  only  needs  g  to  be 
positive  and  to  have  a  suitable  number  of  derivatives  within  a  neighborhood  of  zero. 
Presumably,  one  can  work  out  precise  theorems  with  "stingy  assumptions",  but  this  is  not 
our  intent  here.  It  might  be  interesting  to  investigate  more  severe  departures  from 
normality  than  we  have  considered.  For  instance,  one  should  expect  a  double  exponential 
prior  (g(0)  =  e-l  ^  /2)  to  behave  quite  differently.  The  optimal  boundaries  for  one  severe 
departure  from  normality  is  known:  It  is  easily  seen  that  two  point  prior  distributions  give 
rise  to  sequential  probability  ratio  tests. 

It  is  tempting  to  speculate  whether  Chernoffs  (1965)  theory  for  small  t  (large  s)  can 
be  extended  from  the  setting  of  a  normal  prior  to  the  more  general  setting  of  a  smooth 
prior.  We  do  not  know  whether  this  is  feasible. 
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